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Abstract 

The lepton helicity distributions in the polarized Drell-Yan process at RHIC 
energy are investigated. For the events with relatively low invariant mass 
of lepton pair in which the weak interaction is negligible, only the measure- 
ment of lepton helicity can prove the antisymmetric part of the hadronic 
tensor. Therefore it might be interesting to consider the helicity distri- 
butions of leptons to obtain more information on the structure of nucleon 
from the polarized Drell-Yan process. We estimate the QCD corrections 
at 0(a s ) level to the hadronic tensor including both intermediate 7 and Z 
bosons. We present the numerical analyses for different invariant masses 
and show that the u(u) and d(d) quarks give different and characteristic 
contributions to the lepton helicity distributions. We also estimate the 
lepton helicity asymmetry for the various proton's spin configurations. 



1 INTRODUCTION 

The measurement of the polarized nucleoli structure function gi(x, Q 2 ) by the Euro- 
pean Muon Collaboration pQ in the late '80s has opened the door to the hadron spin 
physics. In the last fifteen years, great progress has been made both theoretically and 
experimentally that has considerably improved our knowledge of the spin structure of 
nucleon. Through these developments, hadron spin physics has grown up as one of 
the most active fields attracting considerable attention. Now our interest has spread 
out to various processes to explore the spin structure of hadrons. In particular, in 
conjunction with new kind of experiments, the RHIC spin project, etc., we are now in 
a position to obtain more information on the structure of hadrons and the dynamics 
of QCD. The spin dependent quantity is, in general, very sensitive to the structure of 
interactions among various particles. Therefore, we will be able to study the detailed 
structure of hadrons based on QCD. We also hope that we can find some clue to new 
physics beyond the standard model through the new experimental data. 

It is now expected that the polarized proton-proton collisions (RHIC-Spin) at BNL 
relativistic heavy-ion collider RHIC [2] will provide sufficient experimental data to 
unveil the structure of nucleon. Therefore it is important and interesting to investigate 
various processes which might be measured in RHIC experiments. One of those will 
be the polarized Drell-Yan process [3]. 

The polarized Drell-Yan process has been studied by many authors both for longitu- 
dinally [UElinilZ] and transversely jHH3IIDllIIllI2] polarized case. The W, Z productions 
from the polarized hadrons are also investigated in this article, we discuss 

the lepton helicity distributions from the polarized Drell-Yan process at the QCD one- 
loop level. The lepton helicity distributions carry more information on the nucleon 
structure ^HIEIEIIE]] than the "inclusive" Drell-Yan observable like the invariant 
mass distribution of lepton pair. Now let us consider, for simplicity, the virtual 7 me- 
diated Drell-Yan process. The subprocess cross section da is written in terms of the 
hadronic and leptonic tensors as, 

da oc (w£ + W A ) (L s » v + L A ^) = W^L S ^ + W A L A » V . 

The anti-symmetric part of hadronic tensor W^ v contains spin information on the an- 
nihilating partons. However, for observables obtained after summing over the helicities 
of lepton or integrating out the lepton distributions, this anti-symmetric part drops 
out. Furthermore, the chiral structure of QED and QCD interactions tells us that only 
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particular helicity states are selected for the q — q annihilation. This observation shows 
that the polarized and unpolarized Drell-Yan processes are governed by essentially the 
same dynamics. On the other hand, if we measure the lepton helicity distributions, we 
can reveal the whole structure of the hadronic tensor. We will show that the u{u) and 
d(d) quarks give characteristic contributions to the lepton helicity distributions. 

The article is organized as follows. In Sec. 2, we reproduce the tree level cross section 
to define our conventions. We present our calculations for the helicity distributions of 
lepton at the QCD one-loop level in Sec. 3. We adopt the massive gluon scheme to 
regularize the infrared and mass singularities avoiding the complexity in the treatment 
of 75. The scheme dependence will be discussed in Sec. 4 and we will change the scheme 
to the MS to perform the numerical studies. We give the numerical results in Sec. 5. 
Finally, Sec. 6 contains the conclusions. The explicit forms for the subprocess cross 
sections are listed in Appendix A. The invariant mass distribution of lepton pair in 
the massive gluon scheme is given in Appendix B which is used to identify the scheme 
changing factor. 



2 POLARIZED DRELL-YAN AT TREE LEVEL 



In this section, we reproduce the tree level result for the helicity distributions of leptons 
in the polarized Drell-Yan process to establish our notation. For the longitudinally 
polarized Drell-Yan process, 



with A (= ±) being the helicity of each particle, we introduce the parton distributions 



where /^(x, +/— ) denotes the distribution of parton type a with positive/negative 
helicity in nucleon A with positive helicity. Based on the factorization theorem, the 
hadronic cross section for the helicity distribution of lepton is given as the convolution 
of the parton distributions with the hard subprocess cross section da ab as, 



N a (Pa, Xa) + N B (P B , X B ) - /(/, A) + 1(1', A') + X , 



tf(x) , Af a A by, 



ffa) = fa& +) + fa& -) , A# (*) = (*, +) - ft (*, -) , 




tfM + X a X A Aff( Xl ) fB( X2 ) + \ b \ B Af b B 



da' 



(X a ,X b ;X) , (1) 
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where 



T 



Q 2 



(Z + Z 



A2 



S (Pa + P> 



The helicity dependent cross section da ab (X a , X^; A) for the subprocess, 

a( Pa , Xa) + 6(p 6 , A 6 ) -> Z(Z, A) + Z~(Z', A') + X , 
is a function of the partonic invariant variables, 

s = (p a + Pbf = (xiPa + x 2 Pb) 2 = XiX 2 S , 

t = {Pa-lf , U = (Pb - If , 



and 



Q* = ( f = {l + l 



l\2 



ZS . 



The tree level polarized cross section for lepton pair production, 



q(p q , Ag) + q(p q , Xg) -> 1(1, A) + Z(Z', A') , 



is given by, 



1 



M 1 



2 • 



2iV c 2 s 

where gZ$2 is the two particle phase space and N c (= 3) is the color factor. The square 
of the tree amplitude reads, 

2 



M 1 



Ana 



where the tree level hadronic and leptonic tensors are defined as, 



W, 



and 



Tr ( uj Xq Pq 7 M Pqlu) , = Tr( ui\ 1 7 M I' j u ) 



IT W%, = 2 [t 2 + u 2 + X q X(u 2 - t 2 )} . 



(2) 



(3) 



In Eq. (jnj) and also below, all momentum p under Tr operator are understood to be p 1 
and uj\ = (l + Xjs)/2. In Eq.(J2J), a is the QED fine structure constant and the quantity 
f XqX depends on the fermion couplings to the photon and Z-boson in the following way, 

Q 2 



r 



e q + Q X « Q X 



9 ^ Q 2 -M 2 + iM z T z 



(4) 
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where e q is the electric charge of the quark in units of the electron charge e, Mz is the 
Z mass, r z is the Z width. The helicity labels X q , A of quark and lepton take ±. The 
quark and lepton couplings to the Z boson are 

q ' Singly COS Ow V ' COs9yj 

where e\ = —1, T 3 is the third component of the isospin and 6w is the Weinberg angle. 
In the center of mass (CM) frame of annihilating quarks, the cross section becomes, 



f x " x (l + cos 2 9 + 2\ q \cos9)5(l- z) , (5) 



da T (X q , \q]\) 7r ( a 

0\,. -As 



dQ 2 d cos 6 ^~^2N C \Q\ 

where 9 is the scattering angle of produced lepton. 

In Eq.(J3J), the third term which depends on the helicities of quark and lepton comes 
from the antisymmetric part of hadronic tensor W^ u . For the observable after taking the 
spin sum of lepton, this antisymmetric part appears in the cross section only through 
the parity violating Z interaction in Eq.(|4"|). Therefore, for the events with "small" 
values of Q 2 , the information from the antisymmetric part will be completely lost. 
Furthermore the chiral structure of Standard Model interaction forces only particular 
helicity states to participate in the process as shown in Eq.(jSJ- These observations tell 
us that for the spin summed over final states and low Q 2 events, the polarized and 
unpolarized Drell-Yan processes are governed by the essentially the same dynamics at 
least for the q q initiated process. 



3 QCD ONE-LOOP CALCULATION 

The principle result of this section will be the polarized Drell-Yan cross section at 
the QCD one-loop level with the helicity of lepton being fixed. At the QCD one-loop 
level, infrared and mass singularities appear and we regularize them by giving a non- 
zero mass k to gluon ^01120] to avoid the complexities from the treatment of 75 and 
phase space integrals. To perform the numerical analyses by using the known MS 
parameterization for the parton densities, we have to change the scheme. However, it 
is well known how to do it [21]. The diagrams to be calculated are given in Fig.l. 
The virtual gluon correction (Figs. la and lb) to this process yields, 

da v da T f a t 

dQ 2 d cos 9 dQ 2 d cos 9 \ it 



In 



2 Q 2 



3 Q 2 7 vr 2 
+ ;t In ^7 - - + 



k 2 2 k 2 4 ' 6 



(6) 
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where a s = g 2 /4:Tc is the strong coupling constant and C F = (iV c 2 — 1)/2JV C for SU(iV c ) 
of color. 

The amplitude for the real gluon emission (Fig.lc and Id), 

q(p q , \) + q{Pg, A,-) -> 1(1, A) + Z(Z', A') + g(k) , 

is given by, 



M 



R 



Q 2 



f X « X U X {l)^Vy{l') 



1 1 



^ Q M A, iPq) e a(^) , 



where e v a {k) is the polarization vector of gluon and T a is the color matrix. We have 
defined q = l + l'. The expression for the square of the amplitude given below has been 
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Figure 1: Parton-level subprocesses contributing to the Drell-Yan process at 0(a s ): 
(a,b) virtual correction, (c,d) real gluon emission (e,f) quark-gluon Compton. 



summed over the spin and colors of unobserved particles. 
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where 



M 



R 



= Sx„-x,N c C F 



(^) \f^wwz. 



w, 



R 



+ 



+ 



i 2 

2 

it 2 
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2(p q ■ /c)Tr(u;A g H^7^) ~ « 2 Tr(cj A 

2(p q ■ tyTriUXsPql,,^) ~ K 2 Tl(uJx q Pql fi Pqlu) 



{2 Pq ■ P q -2p q -k- 2p- q • k) WL + ( Pq , + Pq^W^ 



+ W^ a k a (p qu +Pq U ) 



(7) 



The tree level hadronic and leptonic tensors W^ u , L^ u were denned in the previous 
section and new invariant variables are introduced, 

In Eq.(J7|), we have dropped terms which do not contribute when k 2 — > 0. It should 
be noted that the k 2 terms in the first and second lines can not be neglected since the 
phase space integral produces 1/k 2 singularity |20| . 
The cross section is given by, 

1 



da 



2Nh 



M 



d$3 , 



where <i$3 is the three particle phase space. To integrate out the irrelevant variables 
in Eq.fjHJ), we take the CM frame of qq: 

p, = ^(l, 0,0,1) , Pg = 2 Y (1,0, 0,-1) . 
Parametrizing the momenta of lepton I and lepton pair q by, 



Z M = |1| (1 , sin 6* cosy? , sin#sin<^ , cos#) , 
(f — (<?° ? Q sin 9 cos , q sin 9 sin , q cos ( 



where 



q 



HI 



s + Q 2 - K 2 



|q| 2 = g 02 -g 2 , 



Q 2 



2 q° — |q| (sin 9 sin 9 cos(y3 — 0) + cos 9 cos 9) 
it is easy to write the phase space element d& 3 as, 

|q| Q 2 dQ 2 d cos 9d0d cos 9 dip 



32 (27r) 5 y/s q° — |q| (sin 9 sin 9 cos(0 — 0) + cos 9 cos l 
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After integrating over the angular variables [TO] and dropping terms which vanish when 
K 2 — > 0, we write the inclusive cross section for the polarized lepton in the following 
form: 



da R (X q , Xg; A) _ da RT (X q , Xq] X;z) uR(n ^ „\ , da RF (X q , X q ; A) 



dQ 2 d cos 9 



dQ 2 d cos 9 



F (Q , z) + 



dQ 2 d cos # 



(9) 
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where the first term contains the infrared and mass singularities as well as terms asso- 
ciated with them and, 



da (X q , Xq] A; z) n ( a 

dQ 2 d cos 9 = Xq ^ Xq 2N C [q 2 



2 1 



(10) 



X 



+ 



with 

F R {Q 2 , z) 



(1 + cosfl + z(l -cosfl)) 4 

8Z 2 

(1 -cos# + z(l + cos#)) 4 

'1 2 Q 2 TT 2 



1 + A 9 A)(1 + cos#) 2 + (1 - A,A)(1 - cos0)V} 



{(1 + A,A)(1 + cos0)V + (1 - A 9 A)(1 - cos#) 2 } 



TV 



2ln ^__u (1 _ z) + l±| ln « 2 



K 2 



+ 2(l + z 2 



In (1 — z)\ , 9 , In 2 



The second term in Eq.fjHJ) gives the 0{a s ) finite contribution whose explicit form is 
listed in Appendix A. It should be noted that the first (second) term in Eq.(JTUJ) is just 
the tree level cross section for the q q annihilation with momenta zp q and pq (p q and 
zpq) times z which arises from the difference between the flux normalizations. These 
terms express the processes with collinear gluon emissions. The function F R can be 
understood to be a probability to emit a collinear gluon with momentum (1 — z)p q 
or (1 — z)pq and does not depend on the helicities of quarks because of the helicity 
conservation of QCD interaction. 

By adding the tree level Eq.(jSJ) and virtual Eq.© contributions to Eq.(J§J), the 
double logarithmic infrared singularities cancel out and the cross section from the q q 
initial states becomes, 



da 



T+V+R 



da 



RT 



dQ 2 d cos 9 



+ 



dQ 2 d cos 9 
+ C F ^2(l + z 2 
da RF 



l- 7 --^C F )5(l-z) + ^(pqq(z)\n Q 



4 7T 

/ln(l - z 

1-z 



71 



2(1 + z 2 



\nz 



K 



where 



dQ 2 d cos 9 



'IV 



l + z 2 



+ -5(l-z) 
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which is the DGLAP one-loop splitting functions P qq . To obtain above result, we have 
used the fact, 

da RT (\ q ,\g; X;z) _ da T (X g , Xg-; A) 

dQ 2 dcos6 [ Z) ~ dQ 2 dcos9 ' 
Finally, the contribution from the quark-gluon Compton process (Fig.le and If), 

q(p q , \) + 9% h) - 1(1, X) + 1(1', A') + q(p' q , X' q ) , 

where h is the helicity of gluon, can be calculated in the same way as before. To avoid 
singularities in the physical region, the gluon mass is taken to be — k 2 2(j and the spin 
projection for incoming gluons reads, 



e a (k,h) e* p (k, h) 



—9af3 + i tap-yS 



k~<p\ 

k-Pq 



The amplitude is, 



M 



c 



x u \' q (P q ) 
and its square becomes, 



6 -^j f« X u x (l)Yvy(V) 
1 1 



q + # i> q - 



T a u Xq (p q )e a a (k,h) 



M 



c 



N C C F 



( 47rag s 



jja> W C 



where 




W c 


= i(i 

S 


+ 


f 2u 

\ si 






+ 


H 


with 





u 



~(1 + X q h) I 1 + Tr(uj Xq p' qlll k lv ) 

+ - (i + x q h) + 1(1- x g h) 1 Tr(cu Aq P ; Ta , p g 7v ) 

S b I 



■it 



4 



- j(l - A g /i) 1 + - Tr(^ A9 k^pg^) + j(l + X q h)p p + -(1 - X q h)p qfl p qv 



1 - A g /i 



2Q 2 



) Tr(u Xq (p' q - k) 



7 M P q lu) 



s = (p q + kf , i = (p q - q) 2 , u = (k- q) 2 . 

From the cross section formula for this process, 

1 



da 



c 



2N C (N 2 - 1) 8 
8 



M 



c 



d§? 



we obtain the polarized lepton distribution and write it in the same form as Eq.(JHJ), 

da c (X q , h] A) _ da CT (X q ; A; z) c 2 da CF (X q , h; A) 

dQ 2 dcos8 dQ 2 dcos9 W ,z,A q ,n)+ dQ 2 dcos6 ■ ^ 

The first term, in this case, is given by, 
CT(\ . \, „\ _ / „, \ 2 



da (X q ; X; z) tt I a 
dQ 2 d cos 6 = 2N C I Q 2 " 



f (13) 



cos^(l + cos^ ( (1 + A ^ A)(1 + C0S9)2Z2 + (1 - A ^ A)(1 - COS6)2 } 



and 



F C (Q> ,z,X q ,h) = (^) 



Pg(z; X q , /») ( In ^ + In - 1(1 + M(l - 2*)) 



The quantity is related to the unpolarized and polarized DGLAP splitting functions 
D qg (z) and AP qg ( 



Pqg(z) and AP qg (z) in the following way, 



A 4 , ft) = i (P 9sW - V AP 9sW ) = (1+A ° ft)(1 -f + (1 - A ^ 2 



Equation (fTTTf) again expresses the tree level cross section for the q q annihilation with 
momenta p q and zk where the q is emitted collinearly from the initial gluon. In contrast 
to the qq annihilation case, the function F c does depend on the helicities of initial 
quark and gluon since the emitted antiquark should have the helicity X q = —X q to 
annihilate into the electroweak gauge bosons. The remaining finite contribution da CF 
can be found in Appendix A. 

The cross section for the qg subprocess can be obtained by changing (X q ,h) to 
(—Xq, —h) and \i «-> v in W^ u in the above formulas. 

From these results, the unpolarized cross sections are easily obtained. As a check, 
we have reproduced the results of Ref. [20] by neglecting the Z boson contribution. 



4 FACTORIZATION AND SCHEME CHANGE 

Before proceeding to numerical studies, we must factorize the mass singularities into 
the parton densities in the MS scheme since the two-loop evolution of them which 
should be combined with the one- loop corrections for the hard part, is given in this 
scheme. 
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The renormalized parton densities at the factorization scale \j? F which are relevant 
to this work are written as, 

dy 



q(x,fi F ) 



Aq(x,fj? F ) 



x y 
x y 



x 



q(y)U 1-- + 



2tt 



P 



qq 



x 



ln§ + ^ 



x 



a. 



X 



X 



-HnMr C qg \- 



Aq(y) i l 



x 
V, 



+ 



at. 



2n 



19 



2rr 

fx" 



P 



qq 



In4 + C 



qq 



x 



where C and AC are finite functions depending on the scheme of factorization. Since 
the mass singularity is of universal nature and not affected by the lepton distributions, 
we can find the relation between different regularization schemes by comparing the 
results calculated in respective schemes for the simpler process. We have calculated 
the invariant mass distribution of lepton pair in our scheme and the results are shown in 
Appendix B. The same quantity has been calculated in the dimensional regularization 
scheme in Refs. 0111 • Note that the result depends on also the prescriptions of 75 in 
this scheme. We change our massive gluon scheme to the one adopted in Ref. [Zj which 
can be combined with the known two-loop evolution equations for the parton densities 



in Refs. 



. We find, 



s~1K - 

u qq 



>MS/ 



-c F a + z 2 



In z 



[l - z) 2 + z 2 



1 - z 
In (z(l 



+ 2{l-z) + 



7T 

y 



-7 *(i 



1 



and 



+MS/ 



2z- 1 



]n(z(l-z)) + 



2 y y " 2, 

We must, therefore, factorize these finite terms together with the singular terms from 
the subprocess cross sections. 

The final results for the hard part at 0(a s 
read for qq annihilation from Eq.(|lljl. 



level in the MS factorization scheme 



da 



RT 



dQ 2 d cos 6 dQ 2 d cos 8 



1 



0:.s 



7T 



c h 
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y 



4 )}6(l-z) + 



7T 



P qq (z) ln% 

v H'F 



+ C F {2(l + z 1 



In fl - z) 



-a + z 2 



In z 



+ (1-*) 
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da RF 
dQ 2 d cos 9 ' 

and for qg Compton process from Eq. (jl2j) 
da q9 (i4) da CT a 



(14) 



dQ 2 d cos 9 dQ 2 d cos 9 2n 



pC(z;\ q ,h) <jln^ + ln 



Q 2 , ,„(i--) 2 



z ) 
+ ~(1 - 2X q h(l - z)) 



da CF 

+ dQ 2 dcos9 ' ^ ' 

5 NUMERICAL RESULTS 

We are now ready to predict the helicity distributions of lepton in proton-proton anni- 
hilation. By denoting the scattering angle of lepton in the CM frame of protons with 
momenta Pa and Pb by 9*, the relation between cos 9 in the CM frame of partons with 
momenta x\Pa and x-iPb and cos#* reads, 

x 2 (l + cos6>*) - xAl - cos9*) 

cos 9 = 

x 2 (l + cos6>*) + Xi(l — cos^*) 

The hadronic cross section for the helicity distribution of lepton is finally given by 
inserting Eqs. (|14ll5J) in terms of 9*, into Eq.Q and changing the parton densities into 
the fi 2 F dependent ones. We use the rapidity yi in stead of the scattering angle cos#* 
of lepton in this section, which is defined by, 

1 , 1 + COS0* 

yi = o m 



2 1 - cos 9* 

We fix the total energy S of the proton-proton system to be v^S* = 500 GeV which 
is the planed highest energy at RHIC. Although our formulas can be applied to the 
process for arbitrary Q, we have chosen two values Q = Mz (the Z boson pole) and 
Q = lOGeV as typical examples. For the strong coupling constant, we use the standard 
two-loop form with Aqcd which is set to the value used in the parton densities. The 
following values for other parameters of the standard model are used: 

1 2 n 

a = — — , sin 9w = 0.2315 , 
128 

M z = 91.187 GeV , T z = 2.5 GeV . 



Denoting the helicities of initial protons by Pa,b(±), there are three cases to be analyzed 
corresponding to Pa(+)Pb(—), Pa(+)Pb(+), Pa(—)Pb(—) configurations. 
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We have organized this section as follows. Firstly, we show the numerical results 
for the helicity distributions of lepton with respect to y\. Secondly, we investigate the 
factorization scale and the parton parametrization dependence of the cross section. 



5.1 The Lepton Helicity Distribution 

Here we use the MS parameterization of the parton densities in Refs. j2UI2S] and take 
the factorization scale \ip and the renormalization scale //r to be fip = = Q. The 
scale dependence will be discussed later. 

We give in Fig. 2 the results at Q = M z for both distributions of leptons with 



P A (+) P B (") 



P A (+) P B (+) 



Pa(-) Pb(-) 




Figure 2: The cross sections for leptons with negative and positive helicity in the three 
configurations for the spins of protons with respect to y\ for Q = Mz- 



negative and positive helicities in the three Pa(+)Pb(— ), Pa(+)Pb(+), Pa(—)Pb(—) 
configurations. The unit of the cross section (vertical axis) is [pb/GeV]. The dotted 
line is the total contribution. The solid line is the contribution from the w-quark, 
namely uu and ug(ug) annihilation sub-processes, whereas the dot-dashed line from 
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the o?-quark sector: we write simply d for the d and s quarks. The long dashed (dashed) 
line is the contribution from the tree level uu (dd). We plot in Fig.3 the similar results 
for Q = 10 GeV. Some comments are in order for these results. Firstly, the main effect 
of QCD correction is just an enhancement of the tree level cross section: the K factor 
is 1.5 ~ 1.8. It does not change significantly the shape of the lepton distributions 
from the tree level one. Secondly, we see that the u and d quarks give different and 



P A W Pb(-) P a W P b(+) Pa(-) P b (") 




-2-10 1 2-2-10 1 2-2-10 1 2 

y. 

Figure 3: The cross sections for leptons with negative and positive helicity in the three 
configurations for the spins of protons with respect to y\ for Q = 10 GeV. 

characteristic contributions to the helicity distributions for lepton. In particular, the 
case of the lepton distribution in the Pa(+)Pb(—) configuration for Q = M z deserves 
some explanations. Our numerical results show that in the negative rapidity region, the 
leptons with negative helicities are mainly produced and those come from the w-quark 
annihilation. On the other hand, in the positive rapidity region, leptons with both 
helicities are produced, however the leptons with negative helicities come mainly from 
the ci-quark annihilation subprocesses. These features can be understood intuitively 



13 



by observing the following aspects. (1) The polarized quark distributions in the polar- 
ized proton P(+) roughly satisfy for the relevant scale Q 2 = Mf and the momentum 
fractions of partons X\X 2 = M z /S, 

u{]) » d(i) ~ u(i) » d(t) » u(t I) ~ d(r, I) , 

where j and j denote the quark's spin parallel and anti-parallel to the parent proton's 
spin. This relation implies the dominant subprocesses in the Pa(+)Pb(—) configuration 
to be (i) ua(\)ub, (ii) uaUb{\), (hi) ^a(I)<^b and (iv) Ja^s(I)- (2) From the angular 
momentum conservation, the spin of produced Z boson is aligned to Pa (Pb) direction 
for uaUb and UaUb (^a^b and g^b) annihilations and the negative (positive) helicity 
lepton from the Z decay has higher probability to be produced in the opposite (same) 
direction of Z boson's spin. (3) The third point to be noted is that the V — A coupling 
is larger than the V + A coupling for the quark- Z boson interaction. This suggests 
that among four subprocesses in (1), (ii) and (iii) eventually dominate the process. 
(4) Finally, since the momentum fractions of quarks are bigger than those of anti- 
quarks, the distributions of negative helicity lepton from (ii) and (iii) are Lorentz 
boosted to the negative and positive rapidity regions respectively. The fact above (3) 
also explains that the d (u) quark contribution is larger than that of u (d) quark for 
the Pa(+)Pb(+) (Pa(—)Pb(—)) configuration. For these configurations, the lepton 
distribution is symmetric around yi = as it should be. 

At the lower value of Q (Q — 10 GeV), the Z boson contribution can be safely 
neglected and the electromagnetic interaction becomes dominant. In this case, the 
situation is not so simple as the previous one (Q = M z ). At first sight, only the 
point (3) above should be dropped. If so, the subprocess (i) and (ii) had dominated 
the process which leaded to the similar distributions as those for the Q = M z case. 
However, the distribution takes somehow different shape due to the following reasons. 
For Q = 10 GeV, the partons with very small momentum fractions participate in the 
process (xix 2 = 10 2 /S). Therefore, for the "dominant" processes for the Q = M z case, 
(i) and (ii), the Lorentz boost effect is fairly large and the Jacobian factor J = sin 2 9* 
apparently diminishes the cross section in the yi distributions. Furthermore at this 
scale, all subprocesses seem to contribute equally to the cross section and the partons 
with small x have less information on the spin of the parent proton. One can see, 
nonetheless, the contribution from the antisymmetric part of the hadronic tensor which 
causes the asymmetric distributions of lepton for the Pa{+)Pb(—) configuration. 
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We also estimate the lepton helicity asymmetry A t which is defined by, 

_ da{\ = -1) - da(X = +1) 
1 = da(X = -l) + da(X = +1) ' 

This asymmetry is plotted in Fig. 4 for the three spin configurations for protons. The 
upper three graphs show the asymmetry for the Q = Mz case and the lower three for 
the Q — 10 GeV case. The solid line is the 0(a s ) result and the long-dashed line is the 

P*(+)P B (-) P*(+)P B (+) P*(-)P B (-) 

1 i — i — ' — i — ' — i — ' — i — ' — i — i — i — ' — i — ' — i — ' — i — ' — i — i — i — ' — i — ' — i — ' — i — ' — i — 




Figure 4: The helicity asymmetry in the three configurations for the spins of protons 
with respect to y\. 



tree level one. In these results, we recognize that there are non negligible effects from 
the QCD higher order corrections other than just the enhancement effects of tree level 
results. This asymmetry amounts to around 20 - 30% in the central rapidity region for 
Q = M z . The asymmetry for Q — 10 GeV is much smaller than that for Q = M z . 
This is consistent with the previous observation that partons with small x which do not 
have enough spin information of parent proton, contribute to the process at smaller Q. 
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5.2 The Scale and Parametrization Dependence 

The physical quantities such as the cross sections or asymmetries are obviously inde- 
pendent of the renormalization and the factorization scales. However the truncation of 
the perturbative series at a fixed order induces the dependence on these scales which 
leads to uncertainty in the theoretical predictions. We study the scale dependence of 
the cross section and asymmetry by comparing the results obtained with three different 
choices for the scale: \j? F = 2Q 2 , \xp = Q 2 and \j? F = 0.5 Q 2 . We will also mention the 
dependence on the various parton parametrizations in this subsection. 

We give the results only for the case of Pa(+)Pb(—) configuration with Q = Mz 
since the results for other spin configurations of protons and the value of Q, are found 
to be essentially the same. The scale dependence of the lepton distributions (the upper 




-2-10 1 2-2-10 1 2 

Figure 5: The scale and the parton parametrization dependence. 

two graphs) as well as the asymmetry (the lower left graph) is shown in Fig. 5 where 
the dashed line corresponds to fi 2 F = 2Q 2 , the solid line to fi 2 F = Q 2 and the dotted 
line to ix 2 F = 0.5 Q 2 . In this calculation, we use the parton densities in Ref. |54"| . One 
can see that the variation is less than 10% in the central rapidity region. 
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For the parton parametrization dependence, we calculate the asymmetry with the 
three sets in Ref. [26J since the difference appears more clearly in the asymmetry than in 
the cross sections. The results are shown in Fig. 5 (the lower right graph). The solid 
line, dashed line and dotted line correspond to GS-A, GS-B and GS-C respectively. 
The dependence on the parton parametrizations is similar to the scale dependence 
in size and not so large. The reason will be that the process is dominated by the 
quark and anti-quark annihilation in the RHIC energy region and the gluon initiated 
Compton subprocess, in which the ambiguity of gluon distribution will appear, gives a 
tiny correction to the cross section. 

6 CONCLUSION 

We have presented a complete calculation at the 0(a s ) order in QCD of the lepton 
helicity distributions in the polarized Drell-Yan process. We have numerically ana- 
lyzed the cross section at Q — 10 GeV and Q = Mz and pointed out that the u(u) 
and d(d) quarks give different and characteristic contributions to the lepton helicity 
distributions which deserve some theoretical interests. The QCD corrections mainly 
enhance the tree level cross sections and this fact can explain qualitatively the lepton 
helicity distributions from the various proton's spin configurations. We have defined 
and estimated the lepton helicity asymmetry which amounts to around 20-30%. 

We have also studied the dependence of the cross section and the asymmetry on the 
scale and parton parametrization. Since the qq subprocess is dominant in the RHIC 
energy region, we will not be able, unfortunately, to find clear difference between various 
parton parameterizations which have big ambiguities in the gluon distributions. 

From the experimental point of view, it seems rather difficult to measure the he- 
licity of produced muon and/or electron in the Drell-Yan process. However, if we can 
observe the r lepton produced from the Drell-Yan process and its decay, we will be 
able to compare the experimental data and theoretical prediction presented in this 
paper. We believe that our calculations and numerical studies may provide a comple- 
mentary information to the existing results for the invariant mass and rapidity or 
distributions of lepton pair in the polarized proton-proton collisions. 
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The explicit forms for the finite contributions to the cross sections are listed in this 
appendix. We use the abbreviation c = cos#. 

1. The real gluon emission process. 



APPENDIX A: THE CROSS SECTION FORMULAE 



da RF (X q , Xg] A) 

dQ 2 d cos 6* 




x 



1 - z I [\ r • :( 1 • r))> ~ ~ V Az 

(1 + A g A)(l + c) 2 + (1 - A g A)/(l - cf ( (1 + c + z(l - c)) 2 
+ (l + c + z(l -c)) 4 11 V 4:Z 



) 



) 



+ 




- 2c 2 (3 + Uz + 10z 2 + Uz 3 + 3z 4 ) - c 4 (l - ^) 2 (1 + z 2 )} 



+ 



2z(l-z) 



3(1 + c + z(l - c)) 4 (l -c + z(l + c)) 4 



x {-(1 + z) A {\ + Uz - 18z 2 + Uz 3 + z 4 ) 



+ 2c 2 (l + ^) 2 (1 + Az + 239/ - 56/ + 239/ + 4z 5 + z 6 ) 
+ 8c 4 z(3 - A5z - 151/ - 118/ - 151/ - 45/ + 3/) 
- 2c 6 (l - ^) 2 (1 + 8^ + 71/ + 48/ + 71/ + 8/ + /) 



+ c 8 (l - z)\l -2z + Qz 2 - 2/ + /)} 
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2. The quark gluon Compton process. 



da CF (X q ,h;X) _ vr / c/ 
dQ 2 d cos 6 ~2N C \Q 2 ~ 



V 7T 



X 



'(l-c + z(l + c)) 2 ' 



Az 



A 7 2 pC(7 , n (l + A g A)(l+ C )V+(l-A g A)(l- C )' 

wl ' 9 ' J (l-c + z(l + c)) 4 

x In 

2z(l - *)(! - / - c(l + Qz + z 2 )) 
(l + c + z(l-c))(l-c + z(l + c)Y 

x |((1 - A 9 A)(1 - c) + (1 + A,A)z(l + c)) Pg{z\ X q , h) 



+ 



((1 - X q h)z(l - c) + (1 + A,/i)(l + c + z{\ - c)))| 



+ -(1 + A 9 A)(1 - \ q h)z 2 {1-z- clnz) 
2z(l- z 2 )cPC(z;\ Q ,h) 

+ — — — L,™ . - + g s + \\g a 



where, 
G s = 



x 



G A = 



(1 + c + z(l - c))(l - c + z(l + c)) 

z(l-z) 

12(1 + c + z(l-c))(l-c+ z(l + c)) 4 

"(1 - /) ((1 - \ q h){\ + 5z + 21/ - 7/ - 2/) 

- 2(1 + X q h)(l + 2z- 42/ + 2z 3 + /)) 

- c ((1 - X q h){A + 16z + 61/ + 98/ - 34/ + 2/ - 3/) 

- 4(1 + X q h){l - 3/ - 92/ - 3/ + / 
+ c 2 ((1 - X q h)(5 + 17 z + 65/ + 208/ + 13/ - 17/ - 3/) 

-72(l + A 9 /i)z 2 (l-z 2 )) 

- 2c 3 ((1 - X q h)z{A + 13z + 38/ + 16/ - 2/ + 3/) 

+ 2(1 + X q h)(l - z) 2 (l -2z- 2/ + /)) 

- c 4 (l - /) ((1 - X q h)(5 -7z + 7/ + 13/) 

-2(l + X q h)(l-z) 2 (l-Az + z 2 )) 
+ c 5 (l - X q h)(l - z) 2 (A + 3/ + 4/ + 3/) 

- c 6 (l - A,/i)(l + z){\ - zf{l -z + z 2 ) 

*(!-*) 

4(1 + c + s(l - c))(l - c + z(l + c)) 3 
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X 



(1 - z 2 ){\ + Az + z 2 ) (-(1 - \ q h)z + 2(1 + X q h)(l - z)) 
+ 2c ((1 - \ q h)z(l + 2z + Az 2 + z 4 ) - 2(1 + X g h)(l + 3z 2 - 3z 
+ 2c 2 (l - z) ((1 - X q h)z 2 {l -2z- z 2 ) + (1 + \ q h){\ - /)) 
- 2c 3 ^(l - z) ((1 - \ q h)(l + z + 3z 2 + z 3 ) + 2(1 + X q h)(l - zf) 
+ c\l-\ q h)z(l-z) 2 (l + z 2 ) 



3 z 5 



APPENDIX B: THE INVARIANT MASS DISTRIBUTIONS 

Here we present the invariant mass distribution of lepton pair in the massive gluon 
scheme before factorizing the mass singularity. Comparing these expressions with those 
in the dimensional regularization, we can find the relation between two regularization 
schemes. 

1. The quark and antiquark annihilation process. 



da T+v+R 
dQ 2 



x z 



,-\ q 
1 ■ 



An ( a 



f 



3N C \Q 2 

- 7 ^C F )s(l-z) + ^(p qq (z)]n^, 

An/ 7T \ K' 

'in (l-z)' 



Inz 



+ C F 2(1 + ^ 2 ) —J ^ _ 2 (l + z 2 )f^-2(l-^) 



l-z 



+ 



l-z 



2. The quark gluon Compton process 



3N C \Q 2 



ZX Yn) 



P?MKh) L^ + ln^ 



1 



1 - 2z + 3z 2 + \ q h {3 - 2z - 3^ 2 }) 
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